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Abstract. 1 

Astronomers have discovered many potential black holes in X-ray bi- 
naries and galactic nuclei. These black holes are usually identified by the 
fact that they are too massive to be neutron stars. Until recently, however, 
there was no convincing evidence that the objects identified as black hole 
candidates actually have event horizons. This has changed with extensive 
applications of a class of accretion models for describing the flow of gas onto 
compact objects; for these solutions, called advection-dominated accretion 
flows (ADAFs), the black hole nature of the accreting star, specifically its 
event horizon, plays an important role. We review the evidence that, at 
low luminosities, accreting black holes in both X-ray binaries and galactic 
nuclei contain ADAFs rather than the standard thin accretion disk. 



1. Introduction 

Astronomers have identified potential black holes in a variety of astrophysi- 
cal objects. In binary star systems, consisting of two stars revolving around 
each other under the influence of their mutual gravitational attraction, 
black holes with masses of order several solar masses (M ~ 5 — 20Mq) are 
thought to exist as the dead remnant of the initially more massive of the 
two stars. 

Occasionally, the black hole in such a system accretes (receives) matter 
from the outer layers of its companion and the binary becomes a powerful 
emitter of high energy radiation (hence the name: X-ray binary, XRB) [135]. 
The radiation is ultimately due to the conversion of gravitational potential 
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energy: matter heats up as it falls into the deep potential well of the black 
hole and the hot gas radiates [38]. 

In addition, supermassive black holes (M ~ 10 6 — 10 10 Mq) probably 
exist at the centers of most galaxies. By accreting matter (stars or gas) in 
their vicinity, these black holes may produce the intense emission that we 
observe from quasars and other Active Galactic Nuclei (AGN) [134]. 

Since the discovery of black hole solutions in Einstein's theory of general 
relativity [100], [61], no definitive proof of their existence has been found. 
The main issue is that a black hole remains, by its very nature, concealed 
to the electromagnetic observer and can only be observed through its in- 
direct gravitational influence on the world beyond its event horizon [82]. 
Mass estimates of compact objects in some systems argue for the existence 
of black holes. In addition, several observational signatures have been pro- 
posed as suggestive of a black hole environment. These lines of evidence do 
not, however, constitute definitive proof of the reality of black holes. 

Recently, a class of accretion solutions called advection-dominated ac- 
cretion flows (ADAFs) has been extensively applied to a number of accre- 
tion systems; for these models, most of the gravitational energy released in 
the infalling gas is carried (or advected) with the accretion flow as thermal 
energy, rather than being radiated as in standard solutions [94]. Since a 
large amount of energy falls onto the central object, there is an important 
distinction between compact objects with an event horizon and those with- 
out; for those with an event horizon, the energy is "lost" to the outside 
observer after it falls onto the central object, while for those with a hard 
surface, the energy is re-radiated and ultimately observed [96]. This new 
and potentially powerful method of discriminating between black holes and 
other kinds of compact objects is the subject of this review. 

In §2, we list the best black hole candidates, and briefly review the 
classical techniques and theories used to detect potential black holes. In §3, 
we describe various solutions to the astrophysical problem of gas accreting 
onto a central object, emphasizing the specific characteristics of ADAFs. In 
§4, we describe how ADAF models have been applied to X-ray binaries and 
galactic nuclei and we explain why these studies strengthen the evidence 
for the existence of black holes. 

2. "Classical" Evidence for Black Holes 

2.1. EVIDENCE BASED ON MASS IN X-RAY BINARIES 

It is well known that neutron stars possess a maximum mass above which 
they collapse to a black hole. The precise mass limit is not known due 
to uncertainties in the physics of neutron star interiors; like white dwarfs, 
neutron stars are supported by degeneracy pressure [16], but unlike white 
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dwarfs, nuclear forces play an important role in their structure [53]. The 
equation of state of matter at ultra-nuclear densities is not well known. The 
maximum mass for the stiffest equation of state is ~ 3 M & [58], and the 
limit goes up by no more than 25 % if one includes the effects of rotation 
[39], [123]. Even without detailed knowledge of the equation of state, it is 
possible from first principles to put a limit on the maximum mass. With 
an equation of state satisfying causality (sound speed less than the speed 
of light) and assuming that general relativity is correct, the limit is ~ 3.4 
M Q [116]. Without causality, it reaches 5.2 M [54], [123]. Although exotic 
compact stars with very large masses could theoretically exist [4], [81], it 
is generally accepted by the astronomical community that compact stellar 
objects with masses above ~ 3 M Q are strong black hole candidates. 

In various X-ray binary systems, one can constrain the mass of the 
central object by observations of the spectral lines of the secondary star 
(usually a main-sequence star). The Doppler shifts of these spectral lines 
give an estimate of the radial velocity of the secondary relative to the 
observer. When combined with the orbital period of the binary, this can 
be translated (through Kepler's 3rd law) into a "mass function" for the 
compact object, 

. Mxsin 3 i .... 

s{Mx) = t+w (1) 

The mass function does not directly give the mass Mx of the compact 
star because of its dependence on the unknown inclination i between the 
normal to the binary orbital plane and the line of sight, and the ratio 
of the secondary mass to the compact star mass, q = A^/Mx- However, 
the mass function is a firm lower limit to Mx- Therefore, mass functions 
above 3 Mq suggest the presence of black holes. Additional observational 
data (absence or presence of eclipses, for instance, or information on the 
nature of the secondary star) can help to constrain i or q, so that a likely 
value of Mx can often be determined. The best stellar mass black hole 
candidates currently known are summarized in Table 1. A detailed account 
of dynamical mass estimates can be found in X-ray Binaries [135]. 



2.2. EVIDENCE BASED ON MASS AND VELOCITY PROFILES OF 
GALACTIC CENTERS 

Lynden-Bell's conjecture [70] that AGN are powered by accreting super- 
massive black holes is supported by the intense and energetic emission seen 
in these systems (implying a high efficiency process like accretion onto a 
compact object, cf. 3.1). Observations of superluminal jets in some sources 
(implying a relativistic source) confirm the hypothesis. In addition, rapid 
variability of the luminosity, in particular the X-ray luminosity, is often 
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TABLE 1. The best black hole candidates among XRBs [20], [133]. Type relers to 
Low Mass X-ray Binary (relatively low mass companion star) vs. High Mass X-ray 
Binary. Mass functions f(Mx) and likely masses Mx are given in solar units (Mq). 



Object 


Type 


f(M x ) 


Likely M x 


References 


GRO J0422+32 (XN Per 92) 


LMXB 


1.21 ± 0.06 


> 9 


[7], [15], [36] 


A 0620-00 (XN Mon 75) 


LMXB 


2.91 + 0.08 


4.9 - 10 


[72], [120] 


GRS 1124-683 (XN Mus 91) 


LMXB 


3.01 + 0.15 


5 - 7.5 


[101] 


4U 1543-47 


LMXB 


0.22 ± 0.02 


1.2 - 7.9 


[103] 


GRO J1655-40 (XN Sco 94) 


LMXB 


3.16 + 0.15 


7.02 ± 0.22 


[5], [102] 


H 1705-250 (XN Oph 77) 


LMXB 


4 + 0.8 


4.9 ± 1.3 


[113] 


GS 2000+250 (XN Vul 88) 


LMXB 


4.97 + 0.1 


8.5 + 1.5 


[14] 


GS 2023+338 (V404 Cyg) 


LMXB 


6.08 ± 0.06 


12.3 + 0.3 


[119], [118] 


0538-641 (LMC-X-3) 


HMXB 


2.3 + 0.3 


7-14 


[22] 


1956+350 (Cyg X-l) 


HMXB 


0.24 + 0.01 


7-20 


[46] 



observed. For a time-scale of variability At, causality implies a source size 
less than R ~ cAt, where c is the speed of light. The estimates of R from 
the observed variability are typically small, thus supporting the black hole 
hypothesis [134]. 

In recent years, new lines of evidence have emerged, thanks to improved 
observational techniques. Detailed spectroscopic studies of the central re- 
gions of nearby galaxies provide valuable information on the line-of-sight 
velocities of matter (gas or stars). Since the presence of a massive black 
hole influences the dynamics of matter orbiting around it [10], these stud- 
ies can help discriminate between black hole and non black hole models. 
Current results favor the presence of black holes, but alternative scenarios 
still remain possible [65]. In Table 2 we list the best supermassive black 
hole candidates revealed by these techniques [117], [37]. 

2.3. EVIDENCE BASED ON X-RAY SPECTRAL LINES 

X-ray emission lines are observed in many Seyfert 1 galaxies (a subclass 
of AGN). The lines are probably due to fluorescence by iron atoms in a 
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TABLE 2. Summary of supermassive black hole candidates 
at the centers of nearby galaxies. The first three are the best 
candidates. The likely mass is given in solar units (Mq). 



Galaxy Likely Mass References 



Milky Way (Sgr A*) 2.5 x 10 6 [27], [50], [51], [45] 
NGC 4486 (M87) 3.2 x 10 9 [52], [71] 

NGC 4258 3.6 x 10 7 [83] 



NGC 224 (M31) 


7.5 x 10 7 


[65], [37] 


NGC 221 (M32) 


3 x 10 6 


[8], [126], [127] 


NGC 1068 


10 7 


[47] 


NGC 3115 


2 x 10 9 


[62] 


NGC 3377 


8 x 10 7 


[65] 


NGC 3379 


0.5 - 2 x 10 8 


[42] 


NGC 4261 


4.9 x 10 8 


[35] 


NGC 4342 


3.2 x 10 8 


[37] 


NGC 4374 (M84) 


1.5 x 10 9 


[13] 


NGC 4486B 


6 x 10 8 


[64] 


NGC 4594 


10 9 


[63] 


NGC 4945 


10 6 


[48] 


NGC 6251 


7.5 x 10 s 


[37] 



MCG-6-30-15 ? [125], [33] 



relatively cold gas irradiated by X-ray photons [43], [79]. As matter nears 
a black hole, the velocity of the gas becomes large. Doppler shifts, as well 
as gravitational redshifts, are expected in the observed radiation. In the 
galaxy MCG-6-30-15, a very broad emission line has been observed and is 
believed to originate from gas rotating close to the central mass at a speed 
comparable to the speed of light [125], [33]. It has been claimed that the 
Kerr rotation parameter of the postulated black hole in MCG-6-30-15 can 
be constrained to a value close to unity [23], though this has since been 
challenged [114]. 
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2.4. BETTER EVIDENCE ? 

Due to the nature of a black hole, it is difficult to prove its existence di- 
rectly. Arguably, the best proof would be based on strong field general 
relativistic effects. However, quite generally, these effects can be well imi- 
tated by a neutron star or a more exotic compact object with a radius of a 
few Schwarzschild radii. The exception is the event horizon, which cannot 
be mimicked by any other object, since it is constitutive of black holes. 
Advection-dominated accretion flows may provide unique evidence for this 
intrinsic feature of black holes [96], [91], [88], [92]. 



3. Accretion solutions 

3.1. WHY ACCRETION ? 

Astrophysics deals with a variety of sources that show either persistent or 
intermittent high energy (typically X-ray) emission. Many of these sources 
involve enormous amounts of released energy. Accretion onto a neutron star 
or a black hole is the most efficient known process for converting matter into 
energy (excluding the obvious matter-antimatter annihilation process), far 
more efficient than the fusion of hydrogen into helium [38], [134]. Therefore, 
one is naturally inclined to explain high-energy astrophysical observations 
through the accretion paradigm. 

In these models, the accreting gas radiates (in part or totally) the energy 
that it acquires in falling into the gravitational potential well of the central 
object. For a given accretion rate M (quantity of matter per unit time), 
the gravitational energy available during accretion onto a spherical object 
of radius Rx and mass Mx is 

where G is the gravitational constant. The radiative efficiency of an accre- 
tion flow is defined as its ability to convert the rest mass of particles at 
infinity into radiated energy, 

Lrad /„-. 

where L ra( i is the luminosity radiated by the flow (less than L acc ) and c is 
the speed of light. For a compact object, i.e. an object with a large ratio 
Mx/Rx, we see from Eq. 2 that the efficiency is potentially very high. For 
a neutron star or a black hole, r\ can be ~ 0.1 (in contrast to the ~ 0.007 
efficiency of hydrogen fusion). 
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It is convenient to scale the luminosity and the mass accretion rate in 
terms of the so-called Eddington limit. The Eddington luminosity is that 
luminosity at which the radiation pressure on a spherical flow balances the 
gravitational attraction of the central object: 



Mx 



JEM = 1.3 x 10 38 ( ^ ) erg s" 1 . (4) 



This is the theoretical upper limit to the luminosity of a steady spherical 
accretion flow (but it is usually meaningful for other geometrical configu- 
rations as well). Another useful quantity is the Eddington accretion rate, 
defined here as the accretion rate at which a 0.1 efficiency accretion flow 
emits the Eddington luminosity: 

We use the following notation for the scaled accretion rate: 

m = — . (6) 

M Edd 



3.2. STANDARD SOLUTIONS 

The equations describing the accretion of gas onto a central object are the 
standard conservation laws of mass, momentum and energy, coupled with 
detailed equations for the microscopic radiation processes. These equations 
are highly nonlinear, and allow multiple solutions. One is generally inter- 
ested only in sufficiently stable solutions, since these are the ones that can 
be observed in nature. 

Historically, the first solution to be considered was Bondi spherical ac- 
cretion [11], in which matter without angular momentum accretes radi- 
ally onto the central object. This solution is usually not relevant close to 
the central object (where the bulk of the observed radiation originates), 
since accreting matter in astrophysical settings invariably possesses a finite 
amount of angular momentum and cannot adopt a purely radial accretion 
configuration. 

The second well known accretion solution is the "thin disk" solution in 
which matter revolves around the central object in nearly Keplerian orbits 
as it slowly accretes radially [121], [99], [108]. To reach the central object, 
the accreting matter must get rid of its angular momentum through vis- 
cous transport processes. Associated with the transport of angular momen- 
tum, viscous dissipation heats the accreting gas, and controls, along with 
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the cooling processes, the energetics of the accretion flow. In the thin disk 
solution, the viscously dissipated gravitational energy is radiated very effi- 
ciently. In terms of the efficiency defined in Eq. 3, a thin disk has r] ~ 0.06 
for a Schwarzschild black hole, and n ~ 0.42 for a maximally rotating Kerr 
black hole [26]. (In the case of a neutron star star, r\ ~ 0.2.) Since the matter 
cools efficiently, it adopts a vertically thin configuration where H/R <C 1, 
R being the distance from the central object and H the height of the disk 
as measured from the orbital plane. The thin disk solution has been the 
common paradigm in accretion theory for many years, and has been suc- 
cessfully applied to disks around young stars, disks in binaries and AGN 
disks [38]. 

The thin disk model assumes that, in steady-state (constant m through 
the disk), the local viscous dissipation rate is everywhere exactly balanced 
by the locally emitted flux. The disk is optically thick, which means that 
each annulus of the disk emits a blackbody spectrum (to first approxima- 
tion) . The overall emission from the disk is the superposition of the emis- 
sion from all annuli, which means a superposition of blackbody spectra at 
different temperatures. However, since the emission is dominated by the 
innermost annuli, what emerges is similar to a single temperature black- 
body spectrum. For a detailed review of thin accretion disk theory and its 
applications, see Accretion Power in Astrophysics [38]. 

The thin disk model has serious difficulties explaining the properties 
of some observed systems. In particular, it cannot account for the fact 
that radiation from many systems often spans the entire electromagnetic 
spectrum, with significant emission from radio to gamma rays; this is very 
different from a simple blackbody type spectrum. 

3.3. ADVECTION-DOMINATED ACCRETION FLOWS 

The thin disk paradigm describes the structure of an accretion flow when 
the gas can (locally) radiate all of the viscously generated energy. When the 
gas does not radiate efficiently, the viscously generated energy is stored as 
thermal energy and advected by the flow. In this case, one has an ADAF, 
whose structure and properties are markedly different from a thin disk. 

ADAFs can exist in two regimes. The first is when the gas is optically 
thin and the cooling time exceeds the inflow time (defined as the time-scale 
for matter to reach the central object) [56], [111], [94], [96], [1]. The second 
is when the gas is so optically thick that a typical photon diffusion time 
out of the flow exceeds the inflow time [60], [9], [3]. In both situations, the 
accretion flow is unable to cool efficiently and advects a significant fraction 
of the internally dissipated energy [18]. We focus on optically thin ADAFs, 
since they have seen the most applications to observed systems (see Narayan 
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[87] and Narayan, Mahadevan & Quataert [93] for reviews of ADAFs and 
their applications). 

3.3.1. Optically thin, two-temperature ADAFs: Basic Assumptions 
Important studies of two-temperature accretion flows and advection-dominated 
accretion flows have been carried out in the past [122], [56], [111], but a de- 
tailed and consistent picture has only emerged recently [94], [95], [96], [1], 
offering the opportunity to apply specific models to astrophysical systems. 

A standard assumption in the theory of hot accretion flows is that the 
gas is two temperature, with the ions significantly hotter than the elec- 
trons. In two-temperature ADAFs (2-T ADAFs), the viscous dissipation is 
assumed to preferentially heat the ions. The fraction of the viscous heating 
that goes directly to the electrons is parameterized by a factor 5, usually 
set to 10~ 3 (~ m e /m p ). In addition, the electrons receive energy from the 
ions via Coulomb collisions. When the density is low, however, this process 
is not very efficient. Therefore, only a small fraction of the viscous energy 
reaches the electrons; this energy is usually radiated (though not always 
[84], [75]). The rest of the energy remains in the ions, which are unable to 
cool in an inflow time. The energy in the ions is thus advected with the flow 
and is finally deposited on the central object. If the object is a black hole, 
the energy disappears; if it has a surface, the advected energy is re-radiated. 

The efficiency of the turbulent viscous transport and dissipation in the 
flow are parameterized, as in the standard thin disk theory, by the "a- 
prescription" [121]: the viscosity is written as v = ac^/Qx, where a is the 
viscosity parameter, c s is the sound speed in the plasma and VLk is the local 
Keplerian angular velocity. If we assume that the Balbus-Hawley instability 
[6] is the origin of the MHD turbulence in the flow, we expect the magnetic 
fields and the gas in an ADAF to be in equipartition (i.e. magnetic pressure 
comparable to gas pressure), and the viscosity parameter a to be close to 
0.3 [55]. 

3.3.2. 2-T ADAFs: Properties of the Flow and the Emitted Radiation 
The dynamical properties of ADAFs were first described by analytical self- 
similar solutions [94] ; these were soon replaced by global numerical solutions 
[2], [90], [19], [40], [107], [105]. The flows are geometrically thick (H/R ~ 1) 
and are well approximated as spherical flows with most of the gas properties 
roughly constant on spherical shells [95]. Matter falls onto the central object 
at radial speeds less than, but comparable to, the free fall velocity and 
with angular velocities significantly less than the local Keplerian value. 
ADAFs are radially convectively unstable, but the implications of this for 
the structure and dynamics of the flows are still unclear [96] . 

A self-consistent determination of the thermal structure of an ADAF 
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requires detailed numerical calculations [96], [88]. The poor energetic cou- 
pling between ions and electrons, and the advection of heat by the ions, 
is reflected in their temperature profiles: the ion temperature goes as T; t ~ 
10 12 K/r (r is the distance from the central object in units of the Schwarzschild 
radius) while the electron temperature saturates at T e ~ 10 9-10 K in the 
inner 10 2 — 10 3 Schwarzschild radii. Above a critical accretion rate, m cr j t , 
the density in the plasma becomes sufficiently high that Coulomb colli- 
sions efficiently transfer energy from the ions to the electrons. The gas then 
radiates most of the viscously dissipated energy and is no longer advection- 
dominated. ADAFs therefore only exist for accretion rates below m cr it ~ a 2 
[56], [111], [96], [1]. For m > m cr it, accretion occurs as a cool thin disk. 

Since ADAFs are optically thin and most of the viscously released en- 
ergy is advected, they are significantly underluminous compared to a thin 
accretion disk with the same accretion rate; in fact, the radiative efficiency 
can be as low as r\ ~ 10~ 3 — 10 -4 at low m. The emission by the hot plasma 
comes almost entirely from the electrons. Bremsstrahlung is an important 
emission process. In addition, since the electrons are in the presence of sig- 
nificant magnetic fields and are marginally relativistic, synchrotron emis- 
sion and inverse Compton radiation are also very important. Synchrotron 
emission is the dominant mechanism at radio, infrared and optical wave- 
lengths, with the peak emission occurring at a wavelength that depends on 

l li 

the black hole mass (X pea k M% ) [73]. Comptonization of soft photons 
by hot electrons contributes from infrared to hard X-rays and is a strong 
function of the accretion rate. Bremsstrahlung emission generally domi- 
nates in the X-ray band at low mass accretion rates, but Comptonization 
takes over at higher accretion rates. A significant amount of 7-rays is also 
emitted, resulting from neutral pions (created by proton-proton collisions) 
which decay into very hard photons [74]. Pair creation processes are not 
very important in ADAFs since the radiation energy densities in the flow 
are low [12], [67]. 

3.3.3. 2-T ADAFs : Stability 

The stability of accretion solutions is an important issue in accretion physics, 
since an accretion solution is not viable if it is unstable (by that one gener- 
ally means linearly unstable). For instance, the important 2-T (non advec- 
tive) solution of Shapiro et al. [122] is known to be violently unstable [106] 
and for this reason will probably never be observed in nature. The stability 
of 2-T ADAF solutions has been investigated in the long and short wave- 
length limits [96], [76], [59], [129], [130]. These studies conclude that, since the 
only unstable modes have sufficiently slow growth rates, ADAFs constitute 
a viable solution. 
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4. Applications of ADAFs 

ADAF models have been applied to a number of low luminosity systems. 
They give a satisfying description of the spectral characteristics of the 
source Sgr A* at the center of our Milky Way Galaxy [92], [77], the weak 
AGN NGC 4258 [69] and M87 [115], and of several quiescent black hole 
XRBs [91], [88], [49]. All of these systems are known to experience low effi- 
ciency accretion and the thin disk solution encounters serious difficulties in 
explaining the observed spectral properties. In addition, the ADAF model 
also fits brighter systems with higher efficiencies [86], [30], [29]. 

It is worth mentioning here that previous studies have shown that the 
spectra of some of these systems can be explained by constructing models 
of optically thin accreting plasmas. However, these models are not neces- 
sarily dynamically consistent and do not explicitly attempt to satisfy mass, 
angular momentum and energy conservation. The ADAF solutions solve 
the observational and theoretical problems, for the first time, in a reliable 
and dynamically consistent way. 

4.1. MODELING TECHNIQUES 

Popham & Gammie [107], among others [2], [105], have computed the steady- 
state dynamical structure of ADAFs in the Kerr geometry. Their numerical 
solutions have been used in the spectral models presented here. The dynam- 
ical model provides the radial profiles of quantities such as radial speed, 
angular speed, density, pressure and viscous dissipation. In these models, 
the structure is determined by 3 parameters : a (viscosity parameter), 7 
(adiabatic index of the gas), and / (advection parameter, i.e. the fraction of 
the viscously dissipated energy which is advected). Usually, models assume 
that a constant fraction (1 — (3) = 0.5 of the total pressure comes from 
the magnetic field pressure and that a ~ 0.6(1 — [5) = 0.3. The adiabatic 
index is given by 7 = (8 — 3/3)/ (6 — 3/3) = 1.44 [28]. The parameter / is 
solved self-consistently by calculating the radiation processes in detail and 
feeding back the information to the dynamical solution [30]. So far, studies 
have been restricted to the Schwarzschild geometry, but ultimately, in full 
generality, the Kerr rotation parameter a of the black hole will constitute 
an additional parameter. 

The steady-state energetic structure of the flow is found by solving the 
electron and proton energy equations. The emission spectrum is obtained 
as a byproduct of this procedure, in which synchrotron, bremsstrahlung 
and Compton processes are all taken into account. Since soft photons can 
be scattered off electrons more than once before escaping the flow, the 
Compton problem has to be solved by a global "iterative scattering method" 
for consistency [88]. In the models presented here, the photon transport 
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Figure 1. (a) Spectrum of an ADAF model of Sgr A* (solid line, based on [92]). The 
mass accretion rate inferred from this model is m = 1.3 x 1CP 4 in Eddington units, 
in agreement with an independent observational determination. Dotted lines show the 
spectra of thin accretion disks, at the same accretion rate (upper) and at rh = 1 x 1CP 8 
(lower). Observational data are shown as circles with error bars, and upper limits are 
indicated by arrows. The box indicates the constraints on the flux and the spectral index 
in soft X-rays, (b) Spectra of ADAF models of Sgr A* where the central mass is taken 
to have a surface at 3 -Rschw and the advected energy is assumed to be re-radiated as a 
blackbody. From top to bottom, the three spectra correspond to rh — 1CP 4 , 1CP 6 , 1CP 8 . 
All three models violate the infrared limit. 




Figure 2. Spectrum of an ADAF model of NGC 4258 (solid line), at an accretion rate 
m — 9 x 10~ 3 , with a transition radius Rtrans = 30i?schw Dotted lines show the spectra 
of thin accretion disks, at an accretion rate rh = 4 x 1CP 3 (upper, adjusted to fit the 
infrared points) and rh = 10~ 5 (lower). 
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Figure 3. Spectrum of an ADAF model of A0620-00 (solid line, based on [88]) at an 
accretion rate m = 4 x 1CT 4 , compared with the observational data. The dotted line 
shows the spectrum of a thin accretion disk with an accretion rate rh = 10~ 5 (adjusted 
to fit the optical flux). 




Figure 4- Spectrum of an ADAF model of V404 Cyg (solid line, based on ref. [88]) at 
an accretion rate rh = 2 x 10 -3 , compared with the observational data. The dotted line 
shows the spectrum of a thin accretion disk with rh — 1.8 x 1CF 3 (adjusted to fit the 
optical flux). 
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Figure 5. The broadband simultaneous RXTE (1.3-12 keV) and BATSE (20-600 keV) 
spectra of Cyg X-l observed during the 1996 low-high (upper panel) and high- low (middle 
panel) state transitions. The bottom panel shows a sequence of ADAF models which are 
in good agreement with the observations [29] . 




Figure 6. An ADAF model of J0422+32 (solid line) in the so-called low state (which 
is actually during outburst), compared with the observational data (dots and errorbars) 
[29]. The dashed line shows a thin disk model at the same accretion rate m = 0.1. 
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has been treated as Newtonian, but including gravitational redshift effects. 
Generic ADAF spectra with photon transport in the Kerr metric have been 
computed by Jaroszynski & Kurpiewski [57]. 

The basic ADAF model has one adjustable parameter: the mass ac- 
cretion rate m. The X-ray flux emitted by the ADAF is very sensitive to 
the density and therefore the accretion rate. For this reason, m is usually 
adjusted so that the resulting spectrum fits the available X-ray data. 

The spectral characteristics of some systems are better explained by 
mixed ADAF + thin disk models. The accretion proceeds via an ADAF in 
the inner regions of the accretion flow, and via a thin disk beyond a so- 
called transition radius (Rtrans)- The contribution of the thin disk to the 
overall spectrum is a strongly decreasing function of Rtrans , so that Rtrans 
is a relevant parameter (and therefore is well determined) only in those 
systems where the thin disk contribution to the emitted flux is significant 
(see NGC 4258 and Cyg X-l, for instance). 

ADAF models and the corresponding observational data points and 
constraints are shown in Fig. 1 to Fig. 6 for a number of well-known black 
hole systems. Spectra of thin accretion disks are also shown. While the 
ADAF models are in qualitative agreement with the observations, the thin 
disk models are generally ruled out quite convincingly. 

4.2. SPECTRAL MODELS OF LOW LUMINOSITY SYSTEMS 
4.2.1. Sgr A* 

Sgr A* is an enigmatic radio source at the center of our Galaxy. Dynamical 
mass estimates show that a dark mass M ~ 2.5 x 1O 6 M0 resides in a region 
of less than 0.1 pc in extension around this source [27], [50], [51]. Based on 
this, it is assumed that Sgr A* is a supermassive black hole. Observations of 
gas flows and stellar winds in its vicinity suggest that Sgr A* accretes matter 
at a rate in the range M ~ a few x 10~ 6 - 1O~ 4 M yr" 1 (m ~ 10~ 4 - 10" 2 ) 
[44], [80]. A standard thin disk accreting with this to at a typical efficiency 
i] ~ 0.1 is far too luminous to account for the observed integrated flux. 
Such a model also disagrees with the observed shape of the spectrum and 
is ruled out (by orders of magnitude) by an upper limit on the Sgr A* flux 
in the near infrared (Fig. la). 

The first suggestion that Sgr A* may be advecting a significant amount 
of energy was made by Rees [110], and the first ADAF spectral model of 
this source appeared in Narayan, Yi & Mahadevan [97]. Since then, the 
modeling techniques [85], [84] and observational constraints have improved, 
allowing the construction of a refined model of Sgr A* [92], [77]. The model 
is consistent with the independent estimates of M and M mentioned above, 



16 



K. MENOU, E. QUATAERT AND R. NARAYAN 



and it fits the observed fluxes in the radio, millimeter and X-ray bands and 
upper limits in the sub-millimeter and infrared bands (see Fig. la). 

In evaluating the quality of the fit in Fig. la, it should be kept in mind 
that only one parameter, m, has been adjusted; this parameter has been 
optimized to fit the X-ray flux at 2 keV. The position of the synchrotron 
peak at 10 12 Hz and its amplitude are not fitted but are predicted by the 
model; the agreement with the data is good. There is, however, a problem 
at lower radio frequencies, < 10 10 Hz, where the model is well below the 
observed flux. This is currently unexplained. 

An important feature of the ADAF model of Sgr A* is that the observed 
low luminosity of the source is explained as a natural consequence of the 
advection of energy in the flow (the radiative efficiency is very low, r] ~ 
5 x 1CP 6 ) and the disappearance of this energy through the event horizon 
of the black hole. The model will not work if the central object has a hard 
surface, as demonstrated in Fig. lb. 

4.2.2. NGC 4258 

The mass of the central black hole in the AGN NGC 4258 has been mea- 
sured to be 3.6 x 10 7 M Q [83]. Highlighting the fact that the observed op- 
tical/UV and X-ray luminosities are significantly sub-Eddington (~ 10~ 4 
and ~ 10 -5 respectively), Lasota et al. [69] proposed that accretion in 
NGC 4258 proceeds through an ADAF. They found that the emission 
spectrum and the low luminosity of this system can be explained by an 
ADAF model, provided that most of the viscously dissipated energy in the 
flow is advected into a central black hole. Since then, new infrared mea- 
surements have been made [17] that constrain the transition radius to be 
Rtrans ~ 30 Schwarzschild radii. The ADAF produces the X-ray emission 
while the outer thin disk accounts for the newly observed infrared emission 
(see Fig. 2, unpublished). The refined model is also in agreement with a 
revised upper limit on the radio flux [109]. Maoz & McKee [78] and Kumar 
[66] find, via quite independent arguments, an accretion rate for NGC 4258 
in agreement with the ADAF model. Neufeld & Maloney [98] estimate a 
much lower m (< 10~ 5 ) in order to explain the maser emission in the source, 
but it is hard to explain the observed spectrum with such an m (Fig. 2). 

4.2.3. Other Low Luminosity Galactic Nuclei 

Quasars are luminous high redshift AGN, which are believed to be pow- 
ered by supermassive black holes with masses around 10 8 — 10 9 M [112]. 
Most nearby bright elliptical galaxies are believed to host dead quasars, 
i.e. quasars which have become inactive through evolution [124], [21]. As 
pointed out by Fabian & Canizares [32], however, the nuclei of these el- 
liptical galaxies are much too dim, given the mass accretion rates inferred 
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from independent methods. Fabian & Rees [34] suggested that the problem 
could be resolved if the accretion in these nuclei is occurring through an 
ADAF. This proposition has been confirmed by Mahadevan [73] and mod- 
els have been developed for specific galaxies such as M87 [115] and M60 
[24]. Lasota et al. [69] have similarly argued that low-luminosity LINER 
and Seyfert galaxies have supermassive black holes in their nuclei accreting 
via ADAFs. 

4.2.4. A0620-00 and V404 Cyg 

Soft X-ray Transients (SXTs) are a class of mass transfer X-ray binaries in 
which the accreting star is often a black hole candidate. Episodically, these 
systems enter high luminosity "outburst" phases, but for most of the time 
they remain in a very low luminosity "quiescent" phase. One of the main is- 
sues in modeling black hole SXTs in quiescence is that a thin accretion disk 
cannot explain both the observed low luminosity and the X-ray flux in these 
systems : at an accretion rate low enough to fit the observed luminosity, 
a thin disk will not emit any significant flux in the X-ray band. Moreover, 
the X-ray spectrum will have the wrong shape. Narayan, McClintock & Yi 
[91] argued that the dilemma can be solved by means of the ADAF model 
and demonstrated this by carrying out spectral fits of A0620-00, V404 Cyg 
and Nova Mus 91 in quiescence. Narayan, Barret & McClintock [88] have 
recently refined these models for V404 Cyg and A0620-00 (see Fig. 3 and 
Fig. 4). The models are in good agreement with the observed shapes of 
the X-ray spectra (recall that the X-ray flux level is fitted by adjusting to, 
but the spectral slope is unconstrained in the fit) and they satisfy other 
observational constraints such as optical measurements and extreme ultra- 
violet upper limits. In the optical, the models are somewhat too luminous 
to fit the data; the shapes of the spectra, however, are predicted well (see 
especially Fig. 3, where the position of the peak is reproduced well by the 
model). In contrast, the thin disk model deviates very significantly from 
the data and is easily ruled out. Hameury et al. [49] showed that obser- 
vations of another SXT, GRO J1655-40, are consistent with the presence 
of an ADAF in quiescence. In addition, their model provides a convincing 
explanation for the time delay that was observed between the optical and 
X-ray light curves during a recent outburst [104]. 

4.3. MORE LUMINOUS SYSTEMS 
4.3.1. Spectral States of XRBs 

Narayan [86] suggested that the various spectral states (from quiescence to 
the luminous outburst phase) of black hole binaries (especially SXTs) can 
be understood in terms of a sequence of thin disk + ADAF models, where 
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Figure 7. A comparison between black hole (BH, filled circles) and neutron star (NS, 
open circles) SXT luminosity variations [41] . The ratio of the quiescent luminosity to the 
peak outburst luminosity is systematically smaller for BH systems than for NS systems. 
This confirms the presence of event horizons in BH candidates. 



the thin disk accounts for most of the luminosity in outburst and is grad- 
ually replaced by an ADAF during the transition to quiescence. Recently, 
Esin et al. [30] have developed these ideas with detailed calculations, and 
succeeded in explaining the observed spectral states of the system Nova 
Mus 91 surprisingly well. The same models also appear to explain other 
luminous black hole systems such as Cyg X-l (Fig. 5) and J0422+32 (in 
outburst, Fig. 6) [29]. These studies, for the first time, unify a fair fraction 
of the phenomenology of black hole XRBs. Once again, the existence of a 
central black hole, into which energy is lost, is an essential feature of the 
models. 

4.3.2. Galaxies and Quasars 

A similar unification is also possible in our understanding of the various 
degrees of activity in galactic nuclei. Yi [131] has proposed an explanation 
for the apparent evolution of quasars from a large population of luminous 
objects at redshift > 2 to a much smaller population in the nearby universe. 
He postulates that the mass accretion rate in galactic nuclei has decreased 
over cosmic time. Supermassive black holes in galactic nuclei had m > 
rhcrit at early times and accreted through thin disks with a high radiative 
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efficiency. As rh decreased with time, the accretion flow switched to a low 
efficiency ADAF. This explains the absence of very luminous AGN in the 
local neighborhood of our Galaxy. 

Di Matteo & Fabian [25] proposed that a significant fraction of the hard 
X-ray (> 2 keV) background could be due to emission from a population 
of galaxies undergoing advection-dominated accretion in their nuclei. Yi & 
Boughn [132] discuss a test of the ADAF paradigm in this context; they 
argue that measurements of the radio and X-ray fluxes could provide direct 
estimates of the central black hole masses in many galaxies. 

4.3.3. Neutron Star SXTs vs. Black Hole SXTs 

As mentioned earlier, SXTs are XRBs that exhibit large luminosity vari- 
ations. Most of the time they remain in a quiescent state, in which their 
luminosities are low and the transferred mass is partially stored in a non- 
steady thin disk. Episodically, they experience an outburst event, during 
which their luminosities become very high [68]. During outbursts, which 
correspond to a sudden accretion of the mass stored in the disk, these 
sources reach luminosities of order the Eddington luminosity (cf. Eq. 4). 
From dynamical mass estimates, and observations of "X-ray bursts" (very 
probably thermonuclear explosions on the surface of neutron stars) , one can 
distinguish between neutron star and black hole candidate systems [135]. 
Since black holes are more massive than neutron stars (and the Eddington 
luminosity scales with the mass of the accreting object), black hole SXTs 
are naturally observed in outburst at a higher maximum luminosity than 
neutron star SXTs. 

Assuming that the ADAF model of SXTs in quiescence [91], [88] is 
generic, Narayan, Garcia & McClintock [89] argued that black hole SXTs 
should experience more substantial luminosity changes from quiescence to 
outburst than neutron star SXTs. In quiescence, the central object accretes 
matter through an ADAF : most of the viscously dissipated energy is ad- 
vected by the flow onto the central object. In the neutron star case, this 
advected energy has to be re-radiated from the stellar surface, whereas 
in the black hole case, it is lost through the horizon. Therefore, black hole 
SXTs in quiescence must be much less luminous (relative to their maximum 
luminosity) than neutron star SXTs, which is exactly what is observed in 
nature (see Fig. 7) . This systematic effect confirms the presence of an event 
horizon in black hole SXTs [89], [41]. 

5. Conclusion 

Advection-dominated accretion flow models of black holes in various XRBs 
and AGN are in qualitative agreement with the best current observations. 
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They provide a satisfying explanation for both the spectral characteristics 
and the extremely low luminosities observed in these systems. The models 
require that a very large fraction of the dissipated energy in the flow is 
advected into the central black hole. If a standard star with a surface, and 
not an event horizon, were present at the center of these systems, changes 
in the spectra and luminosities by orders of magnitude are predicted (e.g. 
compare Fig. la and Fig. lb). The success of ADAF models thus constitutes 
strong evidence for the existence of stellar mass and supermassive black 
holes in the Universe. 

We conclude with an interesting observation. If one considers the candi- 
date black hole XRBs listed in Table 1, apart from LMC X-3, all the other 
systems spend the bulk of their time in the ADAF state and only occasion- 
ally switch to the more efficient thin accretion disk state. Similarly, in Table 
2, the three best candidates, Sgr A*, NGC 4258 and M87, appear to have 
ADAFs, and it is plausible that many of the other systems have ADAFs 
as well, except for NGC 1068, MCG-6-30-15, and possibly NGC 4945. If 
ADAFs 3X6} 3S the evidence suggests, very common, there is in principle no 
difficulty finding systems in which to test for the presence of event horizons. 
The challenge at the moment is primarily observational — ADAFs by their 
nature are very dim and one needs instruments with superior sensitivity to 
achieve adequate signal-to-noise. 
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